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The objective of the NASA-sponsored Aerothermal Modeling Program, Phase 
II — Element A, is to develop improved numerical schemes for predicting combus- 
tor flow field. The effort consists of three technical tasks. Tasks 1 and 2 
have been completed. Task 3 is in progress. 

TASK 1 — NUMERICAL METHOD SELECTION 

Task 1 involved the evaluation of various candidate numerical schemes and 
selection of the promising schemes for detailed assessment under Task 2. The 
criteria for evaluation included accuracy, computational efficiency, stability, 
and ease of extension to multidimensions. The candidate schemes were assessed 
against a variety of simple one- and two-dimensional problems. These results 
led to the selection of the following schemes for further evaluation: 

o flux-spline schemes (linear and cubic) 

o Controlled Numerical Diffusion with Internal Feedback (CONDIF) 

To improve the computational efficiency, a direct inversion technique was 
also selected for further testing. In this approach, the continuity and momen- 
tum equations are solved directly, rather than sequentially. 

TASK 2— TECHNIQUE EVALUATION 

Task 2 involved an in-depth evaluation of the numerical schemes selected 
in Task 1. The accuracy was judged by solving test problems for which refer- 
ence solutions are available. The test cases included problems of scalar 
transport, laminar flows, and turbulent flows. These results indicated supe- 
rior performance of the improved schemes. From scalar transport problems it 
was seen that the cubic flux-spline results were more accurate than those from 
the linear flux-spline. However, the cubic spline involved much more computa- 
tional and programming effort and was not considered for fluid flow calcula- 
tions. 

For all the test problems, the linear flux-spline results were more accu- 
rate than the CONDIF results. The flux-spline scheme exhibited mild oscilla- 
tions in the regions of steep gradient. However, it was felt that the presence 
of physical diffusion would tend to diminish these oscillations. 
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To improve the computational efficiency, the flux-spline (linear) was com- 
bined with a direct inversion technique using the Yale Sparse Matrix Package 

(YSMP) [ref 13. Use of such a technique resulted in a factor of 2 to 3 reduc- 
tion in the computational effort compared to the sequential solvers. A summary 
of the Task 2 effort is presented in references 2 and 3. 

TASK 3-3D COMPUTATIONAL EVALUATION 

Task 3, currently in progress, involves the incorporation of the flux- 
spline scheme and direct solution strategy in a computer program for 3D flows. 

Due to the large storage requirement for the LU factorization, it is not 
possible to invert the continuity and momentum equations for the entire 3D 
field. Consequently, a plane-by-plane solution strategy was devised in which 
the cross stream (in-plane) velocities and pressure are solved in a coupled 
manner and the axial velocity is solved decoupled. However, the axial momentum 
and continuity equations are satisfied simultaneously. Such a procedure used 
in conjunction with the power-law scheme [ref 4] for convection-diffusion was 
found to be fast convergent and robust. Work is continuing on the use of the 
flux-spline scheme. 

To demonstrate the accuracy of the flux-spline scheme for 3D flows, results 
are presented for the following two test cases: 

o radial heat conduction in a rotating hollow sphere 
o shear-driven laminar flow in a cubic cavity 

Radial Heat Conduction in a Rotating Hollow Sphere 

The problem is shown schematically in figure 1. A hollow sphere with its 
center located at the origin of a fixed Cartesian coordinate system rotates 
about the x-axis with a constant angular velocity w = w "e x . The radius of the 
the inner surface is rj_, and is maintained at a uniform temperature T^ ; r 2 is 
radius of the outer surface, which is at temperature T 2 . For the case consid- 
ered here, the radius ratio (r 2 /r^) is taken as 2. 

With uniform properties under steady state, the temperature distribution 
is given by: 
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This problem, which is one-dimensional in the radial direction, appears 
three-dimensional if formulated in a Cartesian coordinate system. The calcula- 
tion domain selected is shown as R in figure 1. The calculation domain is as- 
sumed to be fixed in space, so that the material within R has a steady velocity 
field given by: 

V = (w e )X(x e +ye +ze) (2) 
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The exact temperature distribution in Cartesian coordinates is obtained by 
transforming equation (1) to Cartesian coordinates: 
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A uniform 11 x 11 x 11 grid was used to discretize the computational do- 
main. Results were obtained for a range of Peclet numbers (Pe = pwr^/D and 
compared with the power- law scheme. Table I shows the error at the center 
point of the domain. The error has been defined as: 


T - T 

computed exact 

' ’ a - t . > 7 x 100 

max mm exact 


Shear-Driven Laminar Flow in a Cubic Cavity 

The flow situation under consideration is shown in figure 2. Due to symme- 
try considerations, the computational domain extended only half cavity width in 
the lateral (z) direction. 


The flow Reynolds number is 400 and a uniform 22 x 22 x 12 (x, y, z) grid 
is employed for computations. The present results have been compared with the 
solution of Ku et al. [ref 5] obtained using a pseudospectral method (25 x 25 
x 13 mode). This solution has been designated as "REFERENCE" in subsequent 
figures. 

Figure 3 shows the velocity profiles of the u-component on the vertical 
centerline and the v-component on the horizontal centerline of the plane Z = 

0. 5. It is seen that for the same number of grid points the flux-spline solu- 
tion is more accurate than the (lower-order) power- law solution. 

Computations for turbulent flows are in progress. 
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Table I.— ERROR AT THE CENTER POINT OF THE COMPUTATIONAL DOMAIN 


Pe = 1 10 


100 1000 


Error (e) power-law 4.095 x 10 3 
Error (c) flux-spline 1.365 x 10 -3 


1.439 x 10' 2 
4.445 x 10" 3 


2.922 x 10 -1 
1.609 x 10~ 4 


2.756 x 10 _1 
1.177 x 10~ 2 


RADIAL HEAT CONDITIONS IN A ROTATING HOLLOW SPHERE: (a) PROBLEM 
SCHEMATIC; (b) DOMAIN DISCRETIZATION PATTERN (x, = y =2 » r / 3; 

2 2 ^ 3 ) • 
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